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Abstract 



We consider an evolution equation with the regularized fractional derivative of an 
order a € (0, 1) with respect to the time variable, and a uniformly elliptic operator with 
variable coefficients acting in the spatial variables. Such equations describe diffusion on 
inhomogeneous fractals. A fundamental solution of the Cauchy problem is constructed 
and investigated. 
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1 INTRODUCTION 



In this paper we consider equations of the form 

(t,x) -Bu{t,x) = f(t,x), t G (0,T], x e M n , (1.1) 

where O^, < a < 1, is a regularized fractional derivative (the Caputo derivative), that is 



l a) u) (t,x) 



r(i-a) 



dt 
o 



(t - t)" q m(t, x) dr - r a u(0, x) 



;i.2) 



<9 2 " , 



«J=1 



J j = l J 



is a uniformly elliptic second order differential operator with bounded continuous real-valued 
coefficients. 

A strong motivation for investigating such equations comes from physics. Fractional diffu- 
sion equations describe anomalous diffusion on fractals (physical objects of fractional dimension, 
like some amorphous semiconductors or strongly porous materials; see E^] arid references 
therein). In normal diffusion (described by the heat equation or more general parabolic equa- 
tions) the mean square displacement of a diffusive particle behaves like const 4 for t — > oo. 
A typical behavior for anomalous diffusion is const -t a , and this was the reason to invoke the 
equation (1.1), usually with B = A, where this anomalous behavior is an easy mathematical 
fact. For connections to statistical mechanics see also [TT| IT7]. 

It is natural from the physical point of view to consider a usual Cauchy problem, with the 
initial condition 

u(0, x) = u (x). (1.3) 

This setting determines the necessity to use the regularized fractional derivative (1.2). If, 
for example, one considers instead of (1.2) the Riemann-Liouville fractional derivative defined 
similarly, but without subtracting t~ a u(0, x), then the appropriate initial data will be the limit 
value, as t — > 0, of the fractional integral of a solution of the order 1 — a, not the limit value of 
the solution itself (see e.g. fH\ Ho]). 

Note that on a smooth enough function v(t) the regularized fractional derivative can 

1 1 v'(t) 

be written as (D^ a 'v)(t) = J — ^^ r - I n f ne physical literature the expression on 

the right is used as the basic object for formulating fractional diffusion equations. However, in 
order to proceed rigorously, one has either to use (1.2), or to consider, as in [22], an equivalent 
integral equation instead of the Cauchy problem (1.1), (1.3). On the other hand, the above 
expression for on smooth functions is a special case of a generalized fractional differential 
operator of Dzhrbashyan-Nersessyan [B] who studied ordinary differential equations with such 
operators. For example, they showed that the solution of the Cauchy problem D^u — Xv = 0, 
v(Q) = 1, has the form v(t) = E a (Xt a ) where E a is the Mittag-LefHer function. See also |T§| . 

The mathematical theory of fractional diffusion equations has made only its first steps. An 
expression for the fundamental solution of the Cauchy problem (1.1), (1.3) with B = A was 
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found independently by Schneider and Wyss [22] and Kochubei [T3] . It was also shown in [22] 
that the fundamental solution is non-negative, which led later [231 EH El I2Z| to a probabilistic 
interpretation of the equation (1.1). In [22] only initial functions Mo £ 5(M n ) were considered. 
A more general situation was studied in [T3] where u was permitted to be unbounded, with 
minimal smoothness assumptions. There are also some results regarding initial-boundary value 
problems (see |2BII22|). 

For general problems (1.1), (1-3), in ^3] a uniqueness theorem for bounded solutions, and 
an exact uniqueness theorem (for the case n — 1) for solutions with a possible exponential 
growth were proved; see Sect. 2 below. There are also several papers devoted to the Cauchy 
problem for abstract evolution equations (1.1), in which B is a closed operator on a Banach 
space ([121 El El HE], and others). 

An obvious analogy with the classical theory of parabolic partial differential equations [7, 
Uni HE] suggests a wide range of problems for the general equations (1.1) which deserve to be 
investigated. It is natural to begin with the construction and investigation of a fundamental 
solution of the Cauchy problem (1.1), (1.3). That is the aim of this paper. 

More specifically, we will construct and study, under natural assumptions upon the coeffi- 
cients of B, a Green matrix for the problem (1.1), (1-3), that is such a pair 

{Z(t, x; 0, Y(t, x; £)}, (* G (0, T], x, f E W 1 ) 

that for any bounded continuous function uq (locally Holder continuous, if n > 1) and any 
bounded function /, jointly continuous in (t, x) and locally Holder continuous in x, there exists 
a classical solution of the problem (1.1), (1.3) of the form 

t 

u(t, x) = j Z(t, x- 0^(0 d£ + Jd\J Y(t - A, x, ; y)f(X, y) dy. (1.4) 

R" R n 

As in |13| . we call a function u(t,x) a classical solution if: 

(i) u(t,x) is twice continuously differentiable in x for each t > 0; 

(ii) for each ieK™ u(t,x) is continuous in t on [0,T], and its fractional integral 

t 

(t, x) = — i-j J{t - r)- a u(r, x) dr 
o 

is continuously differentiable in t for t > 0. 

(iii) u(t,x) satisfies (1.1) and (1.3). 

The results of this paper (as well as those of the earlier one ^H]) demonstrate a number 
of interesting features of the equation (1.1), which represent a peculiar union of properties 
typical for second order parabolic differential equations (like a kind of the maximum principle; 
see the proofs of Theorems 1,2 in [T3] ) and general parabolic equations and systems. It is 
well known [3 GDI EE] that the latter are characterized by the parabolic weight 2b (b > 1), so 
that the differentiation in t has the same "force" as the differentiation of the order 2b with 
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respect to the spatial variables. In other words, the differentiation in spatial variables has the 
weight 1/26 with respect to d/dt. For the equation (1.1) we have a/2 instead of 1/26, and this 
formal analogy goes through all the constructions and estimates including the description of 
uniqueness and correctness classes. 

In contrast to classical parabolic equations, the fundamental solution Z of the Cauchy 
problem (1.1), (1.3) (even for B = A) has, if n > 1, a singularity not only at t = 0, but also at 
the "diagonal" x = £ with respect to the spatial variables. This causes serious complications for 
the implementation of the classical Levi method [7|nniCSl- The estimates of the iterated kernels 
require a series of additional regularization procedures absent in the classical case. Moreover, 
while the classical heat kernel has a simple expression in elementary functions, in the fractional 
case we have to deal with Fox's H-functions and their asymptotics, which makes even simple 
(in principle) transformations and estimates difficult technical tasks. Note also that, in order to 
consider an inhomogeneous equation (1.1), we have to develop Levi's method twice - separately 
for Z and Y. 

The results can be generalized easily to the situation, in which all the coefficients of B 
except the leading ones depend also on t. Moreover, if we consider only the case of the zero 
initial functions, that is we look only for the function Y, then we can include, without any 
changes, the general case of time-dependent coefficients. The same is true for the case n — 1 
which resembles the classical theory of parabolic equations. In fact, it seems probable that 
the stationarity assumption can be dropped completely. However, then the machinery of Fox's 
H-functions would not be available, and we would have to deal directly with contour integral 
representations of solutions of fractional diffusion equations with coefficients depending only 
on t. Their study is a complicated task in itself, and it looks reasonable to consider first 
the stationary case, in order to identify the main differences between the equations (1.1) and 
classical parabolic equations. 

The main results of this paper are collected in Sect. 2. Sect. 3 contains the information 
regarding H-functions, and miscellaneous lemmas used subsequently. Proofs of parametrix 
estimates are given in Sect. 4 and used in Sect. 5 for substantiating the Levi method for 
the case where n > 2. The case n — 1 is considered in Sect. 6. In Sect. 7 we prove the 
nonnegativity of the functions Z and Y. 



In this section we describe the assumptions on the coefficients and formulate principal results. 
The proofs will be given in subsequent sections. 

2.1. Equations with constant coefficients. Let us begin with the case where the coefficients 
of B are constant, and only the leading terms are present, so that 



We assume that the matrix A = (a^) is positive definite. 

The function Z = Zo(t,x — £) for this case is obtained by a change of variables from the 
fundamental solution found in [22l EI] for B = A. In order to formulate the result, we recall 
the definition of Fox's H-function which will be used systematically in this paper. 



2 MAIN RESULTS 
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Let //, v,p, q be integers satisfying the conditions < z/ < 1 < (i < q. Suppose we have 
also the complex parameters ci, . . . , c p and d\, . . . , d q , and positive real parameters 71, . . . , j p 
and Si,...,S q , such that Pi n P2 = where 

Pi = {s = -(d 3 + k)/S 3 , j = 1, . . . ,//; fc = 0, 1,2, . . .}, 

P 2 = {s = (l- Cj + fc)/ 7i , j = 1, . . . , u; k = 0, 1, 2, . . .}. 
We will need only the case where 
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fe=i 



fe=i 



The H-function 
is defined by the contour integral 



(ci,7i), 
(d 1 ,S 1 ), 



{d^Sfj)] (dft+i, 5fj, + i), (d q ,S q ) 



2tu / E(s)F(s 



(2.1) 



(2.2) 



with 



c(s) = J] r(4 + M, d(s) = JJ r(i - Cfc - lk s) 



k=l 

g 



k=l 



e(s)= J] r(i-4-M, F( s )= J] r(c fc + 7 ^) 

fc=/i+l k=v+l 

(Pi and P 2 are the sets of poles for C(s) and D(s) respectively). The integration contour L is 

an infinite loop running between s = —00 — ia and s = —00 + io where o > max {| ImdA/Sj} 

i- 3: /' 

in such a way that Pi lies to the left of L, and P2 to the right of L. 

Nearly all classical special functions can be represented as H-functions with appropriate 
parameters. The theory of the H-function including its analytic properties, asymptotics, various 
relations, is expounded in [5J 1201 12H] ■ Below we formulate all the results we need. 

Under the above assumptions the function H^(z) is holomorphic in a certain sector con- 
taining the positive real half-axis. 

Now we can write the formula for the fundamental solution Zq\ 



Zo(t,x-g) 



7T 



-n/2 



(det A)V2 



.M=l 



-n/2 



x H 



n 

-r«£^ (y) (*<-fc)(*i-&) 



(l,a) 



(2.3) 
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where A = (A^n is the matrix inverse to (ay). Two different proofs of (2.3) can be found 
(for B = A) in ^3] and [22J (where an equivalent formula is given). In [T3] we used the 
Fourier transform with respect to the spatial variables; the resulting equation is solved using 
the Mittag-Leffler function, and then the inverse Fourier transform is performed on the basis 
of the appropriate formulas for H-functions. The authors of |22] use d the Mellin transform in 
t, the explicit expression of the Green function for the Laplacian, the inverse Mellin transform, 
and integration formulas for H-functions. 

The function Y = Y (t — A, x — y) appearing in the representation (1.4), for this case has 
the form 



Y (t -X,x-y) 



71 



-n/2 



(det A) 1 / 2 



-i,j=l 



-n/2 



x(t- \) a ~ l H w 



12 



1 n 



x j Vii 



[a, a) 

;t,i), 



• (2-4) 



In fact, Yo(t, x) is the Riemann-Liouville derivative of Zo(t, x) in t, of the order I — a (for x ^ 0, 
Zo(t,x) — > as t —>■ 0, so that the Riemann-Liouville derivative coincides in the case with the 
regularized fractional derivative). 

Estimates of the function Z Q and Y , and of their derivatives, are given in the following 
propositions. Denote R = t~ a \x\ 2 . Here and below the letters C, a will denote various positive 
constants. 



Proposition 1 (see [HI), (i) If R>1, 



then 



\D™Z (t,x)\ <Cr^exp (-ar—\x\—) 



\m\ < 3, 



Ja) r. , "("+ 2 > 



a , , 2 



, t } Zo{t,x) <Ct 2 exp y— at a-«|a;|a-a 
(ii) If R<1, x y^O, then 

\D™Z (t,x)\ < Cr a \x\- n+2 -W, \m\ < 3, 
if n > 3, or n = 2, m ^ 0; 

\z (t,x)\ <ct- a [\\o g (r a \x\ 2 )\ + i], 

if n = 2; 



d m Z (t,x) 



dx r 



(2.5) 
(2.6) 

(2.7) 
(2.8) 
(2.9) 



if n = 1. 
(iii) IfR<l,x^0, then 



Z (t,x) 



Ct- 2a \x\- n+2 , ifn>3, 
< <( Cr a [|log(t- a |x| 2 )| + 1] , t/n = 2, 
Cr 3 "/ 2 , i/n=l. 



(2.10) 
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Note that the orders of the singularities at x = in (2.7), (2.8), and (2.10) are precise 
are based on the asymptotic expansions of the H- functions). 

The next proposition contains estimates of the function Y and its derivatives. 

Proposition 2. (i) If R> 1, then 

„ , oc(n-\-m) f I / a , , 2 \ 

\D™Y (t,x)\ < Ct s 1+a exp f-crr^|x|^J , |m| < 3, 

(ii) If R < 1, x ^ 0, n > A, then 

\D™Y (t,x)\ < Cr Q - 1 |a;r n+4 - |m| , |m| < 3, 

(iii) IfR<l,x^0,n = 4, then 

\Y (t,x)\ ^cr^O io g (r a |x| 2 )| + i], 

|£>™y (M)| < CT 2 ^! log(r>| 2 )| + 1], |m| = 2, 

|^y (t,x)| < cr^-^r 1 !! iog(r a |x| 2 )| + 1], |m| = 3. 

(iv) If R < 1, x ^ ; n — 3, then 

\Y {t,x)\ < cr%-\ 

\D x Y {t,x)\ < cr a ~\ 
\D™Y (t,x)\ < cr 0-1 ^!" 1 , H = 2 > 

|D™y (*,a;)| < cr 0-1 ^! -2 , \m\ = 3. 

(v) If R<1, x ^0, n = 2, then 

\Y (t,x)\<Cr\ 
\D x Y (t,x)\<Ct-^-\ 
\D™Y (t,x)\ < Cr a - l [\ log(r a |a:| 2 )| + 1], \m\ = 2, 
\D™Y (t,x)\ < Cr a - 1 \x\- 1 [\ log(r>| 2 )| + 1], |m| = 3. 
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(vi) If R<1, x ^0, n = 1, then 

\D™Y (t,x)\<Ct a \ m = 0,l,2,3. (2.25) 

2.2. The general case. We make the following assumptions on the coefficients of the operator 

B. 

(Bi) The coefficients a^x), bj(x), c(x) are bounded Holder continuous functions on W 1 . 
(B 2 ) The uniform parabolicity condition: there exists such a constant 5 > that for any 

n 

In order to describe regularity properties of solutions of the equation (1.1) with respect to 
the variable t, it is convenient to introduce the class H^ +x [0,T] of such functions (p(t) that 
t^ipit) is Holder continuous on [0, T] with the exponent a + A, A > 0. It was proved in [T3] that 
the problem (1.1), (1.3) cannot have more than one bounded classical solution belonging (for 
each fixed x G W 1 ) to H« +x [0, T], if < A < 1 - a, < fi < A + 1. 

Below we denote by 7 various Holder exponents with respect to spatial variables (without 
restricting generality they will be assumed equal). 

Theorem, a) There exists a Green matrix {Z(t, x; £), Y(t, x; £)} for the problem (1.1), (1-3), 
of the form 

Z(t, x; = Z (t, x - £; + V z {t, x; 0, 
Y(t,x;Z) = Y {t 1 x-bZ) + V Y {t 1 x;t), 

where (Zq,Y ), the Green matrix of the Cauchy problem for the equation obtained by "freezing" 
the coefficients a,ij at the parameter point £ and setting other coefficients equal to zero, satisfies 
the estimates listed in Propositions 1,2, with the constants independent o/£ G W 1 . The functions 
Vz,Vy satisfy the estimates 

(i) Ifn = l, then 

\D™V z (t,x;£)\ < ct(7-™-iW2 f- a r^\x - , 
\D™V Y (t,x;£)\ < Ct a - 1+ ^- m ~ 1)a/2 exp (-af^ \x - 0^) , 

m = 0,1,2. 

(ii) If n = 2, then 

\V z (t,x;0\<Ct^- a exp {-at~^\x - , 
lAcVzft x; 01 < Cr^- Q |x - 0- 1+7 - 70 exp (-ar^\x - 0^) , 

|Vy(t,x;OI < Ct^ _1 exp (-at~ 5=s|x - £|*== J , 
lAJM^OI < C* * +7 ~ 70 exp f-ar— \ x -£\—\ , 
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where 70 < 7 is an arbitrary fixed positive constant (here and below in the formulas 
containing 70 the constant C depends on j ). 

(iii) If n = 3, then 

\D™V z (t,x;£)\ < C^- a |x-£r Hm|+7-70 exp (-at~^\x - , 
\D n ;V Y {t,x-i)\ < Ct^+^-^x - ^-W^-^ 2 exp (-<rt-*5;\x - , 

\m\ = 0, 1. 

(iv) If n = A, then 

\D™V z (t, x- 0| < Ct^^-^x - ^|-2-IH+7-2 7 o exp (- at -2^\ x - , 
\D™V Y (t,x;0\ < Ct^^-^x - ^|-2-^|+7-27o exp (-ar^\x - , 

\m\ = 0, 1; 70 < 7/2. 

(v) Ifn>5, then 

\D™V z (t,x;£)\ < C^- a |x-e| 2 "" Hm|+7 " 70 exp {-at'^x - , 
\D™V Y {t,x;0\ < Ct^+^-^x - £|2-n-M+(7-7o)/2 exp ^_ <Tf - 5 a_j a ._^^ ^ 

|m| = 0, 1. 

(vi) If n>2, \m\ = 2, then 

\D™V z (t,x;£)\ < Ct^- a \x - i\- n+ ^ exp (-o-r^|x - , 
|£>™Vy(t,a;;f)| < C^-^-^r^exp (-af ^ |ar - , 

where 71,72 are some positive constants. 

(vii) T/ie estimates for the fractional derivatives V^, D^-Vy 0/ i/ie functions Z, Y are the 
same as the above estimates for their second order derivatives. 

b) Let u(t,x) be the solution (1-4) of the problem (1.1), (1-3) with a bounded continuous 
function u Q (locally Holder continuous if n > 1) and a bounded jointly continuous and locally 
Holder continuous, in x, function f . If < A < 1 — a, a + A < u, the for any x e M n 

U(;X)EH« +X [0,T]. 

c) The functions Z and Y are nonnegative. 

The rest of this paper is devoted to the proof of the Theorem. As it was mentioned in 
the Introduction, it can be generalized to the nonstationary case, and some lemmas below are 
formulated in a way general enough to cover that more general situation. 
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3 AUXILIARY RESULTS 



3.1. H- functions. We collect here some results regarding H-functions, which will be used below. 
For the proofs and further details see EDI EH] • 

a) Differentiation formulas. 



d T7-20 



(ci,7i) 



-1 rr30 
H 23 



(ci,7i), (0,1) 

(<mo, (4,5 2 ), (i,i) 



(3.1) 



d tt30 



(ci,7i), (02,72) 

(d 2 ,5 2 ), (d 3 , 5 3 ) 



z n M 



(ci,7i), (c 2 ,72), (0,1) 
(rfi,5i), (d 2 ,5 2 ), (d 3 , 7a) (1,1) 



In a similar way 4-H^2 is expressed via etc. 
By (2.2), 



(l,a) 

(f,i), (1,1) 



1 rr(f + s)r(i + s) 



2?rz 



r(l + as) 



2; 



(3.2) 



(3-3) 



We will need this function with z = u(t — r) Q , u; > 0, t > r. We have 

1 /■ r(f + s )r(i + s) _ 



u(t) = H£ 



w(t-r)" 



(1,«) 
(f,l), (1,1) 



2m 



r(l + as) 



c^ s (t-r) as cis. (3.4) 



In order to calculate the Riemann-Liouville fractional derivative (D° + u) (t) 



1 a 



r(i-a) at 



(s) (is, we note that D® + transforms (t — r) as into 



Now we find from (3.4) that 

(D« + u) (t) = {t- t)-*H, 

Similarly 



T(l + as) 
T(l — a + as) 



(t-T) 



u(t-r)- a 



(1 — a, a) 
(f,l), (1,1) 



(see |2H). 



(3.5) 



(t) = (t - r)"- 1 ^ 



uo{t-r)- a 



(a, a) 

(f,l), (1,1) 



(3.6) 



b) Asymptotics at infinity. The asymptotic behavior of H-functions for z — > 00 has been 
thoroughly investigated. We will need only one result of this kind, for a specific class of H- 
functions, for a real argument, with only the leading term of the asymptotic expansion: 



(ci,7i) 
(di,S 



, • • • , \^pi 



7p) 



1), 



(d M , ^ 



const • z {l - a)lp exp 



(3.7) 



where a = ^ c k - J2 d k + ~ P + 1), b = fl 7fe fc fl ^ and P is g iven by (2.1). 

k=l k=l k=l 1=1 
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In particular, from (3.7) we get estimates of some specific H-functions for \z\ > 1: 

< c\z\^=^ exp \-a\z\ ; (3.1 



rr20 
H 12 



(l,a) 

(1,1), (1,1) 



tt20 
H 12 



n+2-2q / 1 

< C\z\ 2< - 2 - a i exp -a\z\ 2 -^ 



(3.9) 



(a, a) 

(f,l), (1,1). 

Similar estimates hold for functions appearing in expressions for derivatives of the above H- 
functions. In all cases the estimates contain the same exponentially decreasing factor while the 
degrees of positive powers of \z\ are different and specific for each case. 

c) Asymptotics near the origin. In order to obtain, for the above class of H-functions, an 
asymptotic expansion near the origin, we write an H-function as 



(ci,7i), 
(di,5t), 



(c P , lp) 



J2 {^Az 



F(s) 



In specific cases, writing Pi explicitly and using well-known properties of the function 
Gamma, we find a required number of terms in asymptotic expansions. In particular, for 
z^ +0 



H 



(l,a) 

(f,l), (1,1) 



z + O (z min ^ > 2 )) , 



r(j-i) 
r(i-o) 

+ n(z 

v 7 ^ -1/2 

r(i-f)' 4 r(i-a) 



if n > 3, 

ifn = 2, (3.10) 
2 ^ + ^&^ /2 + 0(z 2 ) 1 ifn = l 



(here and below we present those terms of the asymptotic expansions which are actually used 
in this work and in [T3]), 



r(f-i) 



H. 



30 
23 



(l,a), (0,1), 
(f,l), (1,1), (1,1) 



-^z + 0(z™^V) 



r(l-a) 

v 7 ^ .1/2 



if n > 3, 



i _VzL r _ 3 v^ -3/2 _j_ n t. 

^ ra-a) rn-^i -r^^ 



2r(i-f) z 1 r(i-a)*- r(i-fi)' 



if n 
if n 



rr40 
-"34 



(l,a), (0,1), (0,1) 
(f,l), (1,1), (1,1), (1,1) 



( r(g-i) 

r(l-a) 



2 + (z min (f > 2 )) 



if n > 3, 
if n = 2, 



r(i-a) 
4r(i-f) z r(i-a)^ 

+ if^) z3/2 + °( ;z2 )' ifn = 1 



(3.11) 



(3.12) 



H 



50 
15 



(l,a), (0,1), (0,1), (0,1) 

(2,1), (1,1), (1,1), (1,1), (1,1) 



( r(f-i) 



r(l-a) 

2loK2 +0(^ 



+ O (^ min (t . 2 )) 



r(l-a) 



^1/2 + 



8T(1-|)~ ' r(l-a) 

ir Jz4_ Z 3/ 2+oM 



if n > 3, 
if n = 2, 



if n 



(3.13) 
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Next we give the expansions for the function appearing in (2.4) and those emerging in the 
course of differentiating (2.4). A different character of their asymptotics is caused by the fact 
that in this case the pole s = 1 appears both in the numerator and denominator of the fraction 
in (2.2) and cancel each other. We have 



tt20 
H 12 



(a, a) 

(f,l), (1,1) 



r(-a) 

2 log 2 , 



2) z 2 + 0(z m ^^), 



r(-f)' 



7i^ 2 + 0(z 3 logz) 

y 3/2 + 20F 



r(-a) 



02 — 

U(f) 



2 log 2 



z 2 + 0(^ 5 / 2 ), 

+ 72 -z 2 + 0(z 3 log z). 
+ 0(^ 2 ), 



,1/2 , 2 v^ -3/2 

^ r(-f) 



if n > 5, 
if n = 4, 
if n = 3, 
if n = 2, 
if n = 1, 



(3.14) 



where 71,72 are certain constants. Next, 



rr 30 
-^23 



(a, a), (0,1) 

(1,1), (1,1) 



2 -^z 2 + o (z™^ 



r(-a) 
3v^ -3/2 



r(-f) 



+ 0(z 3 log z) 
^ 2 + 0(^ 2 ), 



r(-a)' 



r(-a) 



2 72 z 2 + 0(z 3 log z), 



VjL 1/2 _ _VzL 3/2 , Q( 2\ 



2r(f )' 



if n > 5, 
if n = 4, 
if n = 3, 
if n = 2, 
if n = 1, 



(3.15) 



rr40 
^34 



Finally, 



O-50 
-^45 



z 



(a, a), (0,1), (0,1) 
(1,1). (1,1), (1,1), (1,1) 



41X1-2) , 



+ ^mi„(f,3)^ 

4^ + 4 7 i^ 2 + 0(z 3 logz), 



r(-a) 



X-a) 



^(^.)^ ^ 2r( — ) ^ ^ — ^" ^^(^ )? 



Ur(f)' 



(a, a), (0,1), (0,1), (0,1) 
(f,l), (1,1), (1,1), (1,1), (1,1) 



f^f, 2 + 0(^(f.3)), 

-8^-8 7 iz 2 + 0(z 3 logz), 

- < 27y^F 3/2 _ 2 , Q( z 5/2\ 

— \ 4T(-§ )^ rf-a)^ 



if n > 5, 
if ra = 4, 



2r¥|)- 3/2 + ^ 2 + ^ 5/2 )' if ™ = 3 > 



c^-4f^ + 4 72 z 2 + C>(z 3 logz), ifn = 2, 



8T(f) 



r(-a)' 

2 i2££ _ 8l2 z 2 + 0(z 3 log 2), 



2; 



^-m^)Z 3/2 + o(z% 



if n = 1, 

(3.16) 



if n > 5, 
if ra = 4, 
if n = 3, 
if n = 2, 
if n = 1, 



(3.17) 



3.2. Miscellaneous lemmas. Let us fix /? 6 (0, 1). For any x, y G M n , A < t denote 



p(t,x;X,y) 



\x - yj 



1-/3 
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Lemma 1. Suppose that 0<r<A<t<T ; x, y,£G M. n . Then 

p(t, x; A, y) + p(A, y; r, f ) > p(t, x; r, f ). (3.18) 

Proof. Let /(y) = p(£, x; A, y) + p(A, y; r, £). Introducing a new variable y = y — x we can 
write 



/(y) = /i(y) = ( ^apJ +v (a-^< 



Denoting 



f \ i-/3 /| r _ |£ _ x | 



MO = ( J +( m V ) > ^>°' 



1-/3 



1-/3 



we find that 

min f(y) = min f\(y) > min/ 2 (r). (3.19) 



r>0 



The derivative /^(O i s positive for r > |£ — x\; thus /2( r ) attains its minimum value on the 
interval [0, \x — £|]. A standard search for its possible local extremum points yields a single 



point, in which the function takes the value ( -j- p ) . On the other hand, 



1-/3 



M» " (tt-^V"'. Ml*-*!)" ' K 



1-/3 



l£ — ^1 

whence /^(r) > ( -; -» ) . Together with (3.19) this implies (3.18). □ 



(t-< 



The next three lemmas establish integral inequalities involving the function p. We keep the 
notation of Lemma 1. 

Lemma 2. For any \x > 0, e G (0, 1) 

[(t - A)(A - r)]~^ | exp{-p[p(t, x; A, y) + p(A, y; r, 0]} 

< C(e)(t - r)-^ exp{-(l - £ )pp(t, x; r, 0}- (3.20) 

This lemma is a consequence of Lemma 1. See the proof of Lemma 5.1 in |J| where a similar 
result was obtained for n — 1. 

Lemma 3. For any ai, ct2 > 0, any &i, 6 2 > 0, suc/j i/iai &i + 62 < ^ ; an d any p > ; e G (0, 1) 

t 

y (t - A)" 1 " 1 ^ - r)^- 1 y |x - y| 6l -"|y - £| b2 ~" exp{-p[p(t, x; A, y) + p(A, y; r, 0]} dy dA 

T R™ 

< C(e, 01, a 2 , 61, 6a)(* - r) 91 ^- 1 ^ - ^ +b ^ n exp{-(l - e)pp(t, x; r, 0}- (3.21) 
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Proof. Writing 





4 

— LJ *k Wlieie 

k=l 












Vi 


= {y e W 1 : \x 


-y\>(t- 




\y 


-ei>(A- 




v 2 


= {y Gl": \x 


-y\>(t- 




\y 


-ei<(A- 


r« 




= {y Gl": \x 


-y\<(t- 




\y 


-ei>(A- 


rf} 


Vl 


= {y Gf": \x 


-y\<(t- 


Af, 


\y 


-ei<(A- 





we decompose the left-hand side in (3.21) into the sum of four integrals Jj, J 2 , J 3 , J 4 . 
If y e 14, then |x - y| bl - n < (t - A)^ 1 "™), |y - £| &2 ~ n < (A - r)^ 2 ""), so that 



t 



J x < J(t- X) ai+(3bl -\\ - r) a2+l3b2 - 1 [{t - A)(A - r)]- n/3 

T 

x J exp{-fi[p(t,x;X,y) + p(X,y;r,^)]} dydX, 

and by Lemma 2 



Ji<C(t-rr +ffl2+ ^ +b -^- 1 exp{-(l-|) / ip(t,a;;r,0}. (3.22) 



Since 



(t _ r) (6 1 +6a-^ exp{ _£ //p(t)j; . T)e)} 

= [p(t, x- r, tff-Mn-h-bz) e xp{-| M p(t, s; T , £)}\x- ^ +b >- n <C\x- C\ b ^~ n , 

the inequality (3.22) means that J\ does not exceed the right-hand side of (3.21). 

In order to estimate J 2 , we consider two distinct cases. Suppose first that \x — £| > 2(t — r) 13 . 
Then \x — y\ > (t — r) 13 , so that by Lemma 1 

t 

■h < (t-T)^- n Uxp{-fxp(t,x;r^)} J (t-Xr^iX-rr-UX J \y - dy 

|2/-?I<(A-t)/3 

= C(t - r)° 1+aa_1 (t - T )^ bl+b2 -^ exp{-/ip(t, x- r, f )} 

< C(e)(t - r)* 4 *" 1 ^ - it +b2 - n exp{-(l - x; r, 0}- (3-23) 

If |z - f | < 2(t - r)' 3 , then \x - y\ < 3(t - rf, 

i 

J 2 < exp{-fip(t, x- r, £)}/(* - A)"- 1 (A - rr- 1 dX J \x - y\ b ^ n \y - i\ b ^ n dy, 

\y-Z\<3(t-T)0 

and by Lemma 2 from Chapter 1 of jTU] we get 

■h < C(t - r) ai+a2_1 |x - £\ bl+b2 ~ n exp{-p,p(t, x; t,£)}. 
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Together with (3.23) this yields the required bound for J 2 . 
The estimate for J3 is proved similarly to that of J 2 . 

Let us consider J4. It follows from the definition of V4 that J4 = if \x — £| > 2(t — r) 13 . 
Therefore we may assume that \x — £| < 2(t — ry and note that 

V,c{yeR n : \y-i\<2{t-rf}. 

Then 

t 

J4 < exp{-pp(t, x- T,£)}f(t- \r-\X - r)^- 1 dX J \x - y\ b ^ n \y - ^ n dy. 

\v-e\<2{t-r)f> 

Using again Lemma 2 from Chapter 1 of jTU] we find that 

J 4 < C(t - r) ai+ai, - 1 \x - ^\ b ^~ n exp{-fi,p(t,x;T,€)}, 

as desired. □ 

Remark 1. If 61 + 62 > n, then the assertion of Lemma 3 still holds if the factor \x — £\ bl+b2 ~ n 
is omitted from the right-hand side of (3.21). The proof is similar; we have only to use the 
appropriate estimate from the same lemma of [TO] . 

Lemma 4. I/O < C2 < c\ < 1 ; 7 > 0, then 

J \x - ?/r n+7 exp{-/i[cip(t,x; X,y) + c 2 p(X,y;r,£)]} dy 

< Cexp{-ne 2 p(t,x;T,£)}(t- Af 7 (3.24) 

where the constant depends on c 2 — c\ . 

Proof. Let us write the integral in the left-hand side of (3.24) as the sum I\ + I 2 of two 
integrals corresponding to the decomposition W 1 = W\ U W 2 where W\ = {y G M. n : \y — x\ < 
(t - X) 13 }, W 2 = {y E W n : \y-x\> (t- X) p }. 

Using Lemma 1, we get 

h < exp{-pc 2 p(t,x; r,f)} J \x - y\~ n+1 dy = C(t - A)^ 7 exp{-pc 2 p(t, x; t, £)}. 

\y-x\<(t-\)0 

Next, again by Lemma 1, 



I 2 = J \x-y\ n+1 exp{-p[cip(t,x;X,y) + C2p(X,y;T,C))}dy 

\y-x\>{t-X)P 

< (t - Xy n(3+ ~<P exp{-pc 2 p(t,x;r,£)} J exp{— ^(ci - c 2 )p(t, x; A, y) dy] 



<{t- Af 7 exp{-pc 2 p(t, x- t, f )}, 
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as desired. □ 

The next two lemmas give estimates of the iterated kernels 

t 

K m (t,x\T,£) = J dX J K 1 (t,x;\,y)K m _ 1 (\,y,T,£,)dy, m > 2, (3.25) 

T R n 

where K 1 = K is a given kernel. We will treat the cases n > 2 and n — 1 separately. 
Lemma 5. Let n > 2. Suppose that K(t,x;r,£) is a continuous function on 
V n = {(t,x;r,0 : x, f E R n , x ^ f, < r < t < T} , 

| x; r, 1 <A 1 (t- rf^x - ^ n+ ^ exp{-^(t, x; r, 0} (3.26) 
where A 1 ,/i > 0, < u ,Ui,P < 1, and nzz-f 1 no£ an integer. Then the series lZ(t,x;T,£) = 

oo 

K m (t,x;T,^) is absolutely and uniformly convergent on V n H {(i, x;r, : t — r > 5 > 

m=l 

0, \x - £| > 5 > 0}, and 

\K(t, x; r, 1 < C(t - rf^x - exp{-//p(t, x; r, £)}, (3.27) 

with some p* G (0,/i). 

Proof The bounds for the kernels (3.25) are obtained by two stages. First we use Lemma 
3. We find, for a small e > 0, that 

\K 2 (t,x; r,0l < C 2 (e)(f - rf^x - fl-* 2 " 1 exp{-p(l - s)p(t,x; r,0}, 
and so on, so that 

|tf s (t, x; r, 1 < C fl (e)(t - r)'^" 1 ^ - CP"^ 1 exp{-/x(l - (s - l)e)p(t, x; r, 0}, (3.28) 

whenever — n + sz/i < 0. Let s* = min{s : —n + sv\ > 0}. Then by Remark 1 

\K S * (t, x; r, 0| < C (f " ^) a(,_1 exp{-/i*p(t, x; r, 0}, (3.29) 

where a = s*u /3, p* = [1 — (s* — 

For all the next iterations, we use (3.26), (3.29), and Lemma 4. Using the identity 



J(t- \) ai -\\ - rf 2 - 1 d\ = B(a u a 2 )(t - r) 

T 

we get by induction that 
\K s * +m (t, x- r, 01 < C Q C? \ J] B (("o + vx)P, a + (is + Vl ){j - 1)) 



Q1+02— 1 



x (i - r ) m (-o+-i)/3+«o ex p{-//p(£, x; r, 0}- (3.30) 
17 



r(ao)fr((z/o + z/ i)/3)l m_1 

The product in the right-hand side of (3.30) equals — — -— -. This implies 

T(a + (m - l)(z/ + vi)f3) 

the convergence of the series for 1Z. The estimate (3.27) follows from (3.28) and (3.30). □ 

Let us consider the iterations (3.25) for n — 1, with the estimates coming from the Levi 
method for this case. 

Lemma 6. Let n = 1. Suppose that the function K(t,x; r, £) is continuous on 
V 1 = {(t,x;r,Z) : x,^ E R 1 , x ^ ^0 < r < t < T} , 

and satisfies the inequality 

\K(t, x- r, 01 < A(t - rf^- 1 exp{-/ip(t, x; r, £)} 

oo 

where Ai, p > 0, < 7, /3 < l.Then the series TZ(t,x;r,^) = ^ K m (t,x;r,^) converges 

771=1 

absolutely and uniformly on V\ D {(t, x; r, £) : t — r > 5 > 0} ; and 

\K(t, x- r, 1 < C(t - r)^" 1 exp{-//p(t, x; r, £)}, (3.31) 
wiz;/i some p* E (0,/i). 

Proof. Using Lemma 2, we find for any e > that 
t 00 
|tf 2 (f,x;T,0l<^? jKt-X^X-r^-^dX J exp{p[p(t,x; X,y) + p(X,y;r,0}dy 

T — OO 

< A 2 (e) (f - r) 2 ^' 3 - 1 exp{-Ml - e)p(*, s; r, 0}- 
Repeating the procedure we obtain the inequality 

|#m(t, ^; r,Z)\< A m {e) (t - tT^- 1 exp{-//(l - (m - l)e)p(t, x; r, 0}- 
'1 + 13' 



Let m* = — — + 1. Then 
. Pi . 

\K m *(t,x;r,$,)\ < An'(e)exp{-^(1 - m*e)p(t,x;T,£)}. (3.32) 
Setting e = , p > 2, p* = /i(l — p" 1 ), A = max A m (e), we can rewrite (3.32) as 

pm* l<m<m* 

|K m *(i,x;r,0l < A exp{-p*p(t,x;r,£)}. (3.33) 

Now, for m > m*, we proceed in a different way in order to preserve the exponential factor 
p* in the estimates of further iterated kernels. By Lemma 1, 



1 00 

\K m , +1 (t,x-T,0\<A 2 Jit-XyP-P-'dX J ex P {p[p(t,x;X,y) + (l-p~ 1 )p(X,y;r^)} 



dy 



T 

t 



< AlF exp{-p*p(t, x- T,£)}J(t- X^- 1 dX = AlFB^/3, l)(t - r)^ exp{-p*p(t, x; r, £)}, 
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00 f f J ' 1 1 
where F = J exp ^ M w3 j 

Similarly, we prove by induction that 
\K m * +j (t,x; r,0\< A (A F) j (t - r)^ j \l[B(^, 1 + S1 (3) \ exp{-//p(t, x; r,0}- 



s=0 



Using the identity 

f = o r(l + 7/3j) 

we see that the series for 7?. is majorized by the convergent series 



This implies (3.31). □ 



4 PARAMETRIX 

4.1. TTie function Z . The parametrix kernel Z (t,x — £,0 is defined by the formula (2.3) 
where (A^) is the matrix inverse to the matrix (f%(C)) of the leading coefficients "frozen" at 
the parametric point ( e 1". The estimates of Z proved in ^3] for the case of the constant co- 
efficients ciij and collected in Proposition 1, remain valid for the parametrix, with the constants 
independent on (. 

Since properties of Zq are different for n > 2 and n = 1, it is convenient to treat these cases 
separately. Therefore we assume in this and the next sections that n > 2. The case n = 1 will 
be considered in Sect. 6. 

The behavior of Z for R > 1 and R < 1 is described in Proposition 1 separately. However 
it is possible to write equivalent unified estimates: 

\D™Z (t,x-£;0\ < Cr Q |x-^|-™ +2 -Hexp{-ap(t,x;0,O}, H < 3, (4.1) 



i a) Z (t,x-t;() < Ct~ a \x - i\~ n exp{-ap(t,x;0,O}, H < 3, 



(4.2) 



where 



p(t,x; r,0 



2 

2-a 



,(t-T)«/ 2 ( 

(this corresponds to the notation of Sect. 3.2 with /? = a/2). If n = 2, then (4.1) is valid only 
for \m\ ^ 0. For the opposite case 

|^o(t, x - £; C)| < CT" lo S ( r > - ^) I + !] exp{-ap(t, x; 0, 0}- (4.3) 

The constants C, a > in (4.1)-(4.3) do not depend on (. Of course, the constants may be 
different from the ones in (2.5)-(2.10). 

The next proposition gives similar estimates for differences of values of Z and its derivatives 
corresponding to different values of the parameter (. 
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Proposition 3. For any y, (" G W l , < t < T, 

\D™Z (t, y ;(')-D™Z (t, y ;(")\ 

< Ct' a \C - ("|7| y |-n+2-|m| exp {_ a (r a / 2 |y|) ^ }, 

Proof. Let m = 0. Denote 

n 

By our assumptions 

|2l(y,0-2t(y,C")l <C|C'-CTM 2 , 



|m| < 2. (4.4) 



[det A(C')] V2 - [det A(C")}^ <C\('- CV- 

According to (2.3), we have to use the estimate (3.8) for the function Hf®, and also to find 
an estimate for a difference of such functions. Let 



//xi 1/2 



<p(s) = s- n/2 H- 



LdS 



By (3.1), 



n 



iff {s) = ~^- l H 



1 cj20 
12 



uos 



(l,a) 



(1,«), (0,1) 
(f,l), (1,1), (1,1) 



For large s, (3.7) and (3.8) give the estimate 

n+2 n -i ra+2 / 1 \ 

Iv 3 ( s )l < Cs 2 ' 2 -") 2 cj 2 < 2 - q ) exp [ -o-|ws| 2 — J , 
which implies the bound 

\<p (2l(y,0) " V »,C"))I < Cr¥|^ - CV exp (-(rt-^|y|^=) . 
For small s, we have 

|<//(s)| < Cus'* 

by (3.10) and (3.11); note that the logarithmic terms emerging for n = 2 are cancelled. 



As a result, we find that 
\Z (t,x-Z;C)-Z (t,x-Z,(")\< 



'ClC'-CT^exp (-<7r^|j/|^=) if r Q |y| > 1, 

.c|C-C"l 7 i- a |yr n+2 , ifr Q | y |<i, 



which is equivalent to the estimate (4.4) with m — 0. 

For the first and second derivatives the proof is similar, though somewhat cumbersome - 
one has to use the asymptotics (3. 7)- (3. 12) for the H-functions; again for n = 2 the logarithmic 
terms are cancelled. □ 
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4.2. The function Yq. The estimates for the function Yq are given in Proposition 2 for the 
case of constant coefficients. They carry over to the kernel Y (t — X,x — £; C) defined by (2.4) 
with the coefficients depending on ( as above. The proof is a direct consequence of the 
formula (2.4) and properties of H-functions, in particular the differentiation formulas (3.1), (3.2) 
etc, the asymptotic relations (3.7) and (3.14)-(3.17). The calculations are simple but tedious, 
especially for n < 4, since we have to take into account several terms of the asymptotics of 
H-functions near the origin. 

Just as for the function Z , it is desirable to obtain unified estimates of Y and its derivatives 
which are valid for all values of the independent variables. Here the problem is a little more 
complicated because the behavior of Yq is different for different values of the dimension n. 

In particular, for n = 2 we have 

\D™Y (t - T ,x-t,Q I < C(t - T )-f 0+M)+«-i exp{-ap(t, x; r, £)}, H = 0, 1; (4.5) 



\D™Y (t-T,x-t;0\ 

< C(i-r)- a - 1 |a;-f| 2 - |m| 



log 



k-CP 



(t-T) 



1 



exp{-ap(t,x;r,£)}, |m|=2,3. (4.6) 



These estimates follow immediately from (2.11), (2.21)-(2.24). 
For n — 3, 



\D™Y (t - r, x - £ C) I < C(t - T )-f («+M)+«-i cxp {- a p(t, x; r, £)}, H = 0, 1, (4.7) 



and 



\D?Y (t -t,x-£;C)\ < C(t- r)— 1 \x - ^"M exp{-ap(t, x; r, £)}, H = 2, 3. (4.8) 
Similarly, for n — 4, 

|£>™y (*-T,x-e;C)l 



< C^-r)-"-^- 1 



log 



k-CP 



+ i 



exp{-crp(t,x;r,0}, H =0,1,2, (4.9) 



|^y (i-r,x-e;C)l 

< C{t-r)- 2a - l \x-i\- 1 



log 



k-cp 



exp{-(7p(i,x;r,0}, H = 3. (4.10) 



(f - r)« 
Finally, if n > 4, then 

|^y (< - r,x - C; C)l < C(t - r)- a - l \x - Z\- n -W +4 exp{-ap(t,x; r,C)}, \m\ < 3. (4.11) 

Just as in Proposition 3, the estimates for the differences D™Y§{t — r, x — £; CO — -D™^o(^ ~~ 
r, x — £;C") have the following form: the right-hand sides of (4.5)-(4.11) are multiplied by 
IC — C"l 7 - F° r example, if n > 4, then 

\D?Y (t -r,x-i; C) - £>™y (* t, x Ci CO I 

< C(t-r)- a - 1 |C / -C"rk-^r n ~ H+4 exp{-<rp(t,x;r,0}, M < 3. (4.12) 
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4.3. Integral identities. It follows from the construction of the function Zq that 



J Z Q {t,x-S;Qd£=\ 
(see 1.13]). Next, from (2.4) we see that 



(4.13) 



j Y {t-r,x-ii-Qdt={t-TY- l Tr- n l 2 J \y\- n H 



n u20 
12 



z (t-r)- a \y\ 2 



C(t - r)"- 1 J r- x Ef 2 



C x {t-r) 1 



a-l / o -lff20 
6 -"12 



-(t-T)~ a r 2 
4 V ; 



(a, a) 

(1,1) 

(a, a) 

(f,l), (1,1). 



(a, a) 

(f,l), (1,1) 



-Ci^-r)"- 1 lim H. 



2:-! 



-(t-r)- a s 



(a, a), (1,1) 

(0,1), (f,i), (1,1) 



so that 



We have used the integration formula 



(4.14) 



1 tj20 

12 



as 



[a, a) 

;t,i), (1,1) 
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at 



(1,1), (a, a) 
(f,l), (1,1), (0,1) 



(pOj. formula 1.16.4.1), the identity connecting and iff^ ([20 , 8.3.2.9), and the asymptotics 
(3-7) of tff 3 . 

4.4. Further estimates. Below we will need also the first time derivative of the function 
Z (t,x;(). We find from (2.3) and (3.1) that 



dZ (t,x;C) an-^H- 1 



dt 



(detA)V2 



-n/2 



(1,«), (0,1) 

(1,1), (1,1), (1,1) 



Using the asymptotics (3.7) and (3.11) of the function we can obtain an estimate of 
For example, if n > 3, then 



dZ 
dt ' 



dZ (t,x;C) 



dt 



< ^Ct- 1 -? exp {-a(t-«ft\x\)&} , if t— /2 |x| > 1, 

cr a - l \x\- n+2 , iit- a ' 2 \x\<i. 
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We can also use a roughened unified estimate 
dZ (t,x;C) 



dt 



< Cr l \x\- n e*v{-a(t- a l 2 \x\)^y (4.15') 

dZ (t,x;C) dZ (t,x;(") 



whose 



As before, we get also an estimate for the difference 

dt dt 

upper bound is the expression in the right-hand side of (4.15) or (4.15') multiplied by |C'~C"I 7 - 
Since 

dZ (t,x-S;C') 



I 



dt 



by virtue of (4.13), the above estimates imply also the estimate 



dZ (t,x-Z;S) 
dt 



< Ct~ 1+ ^ 



(4.16) 



dZ^it x~ x) 

obtained by subtracting the (zero) integral of ^— - — and using the estimate for the dif- 



ference of the derivatives. 



dt 



5 THE LEVI METHOD (n > 2) 

5.1. The scheme. We look for the functions Y, Z appearing in (1.4) assuming the following 
integral representations: 

t 

Z(t, x- = Z (t, x-£;0 + JdxJ Y (t - A, x - y; y)Q(X, y; dy; (5.1) 



Y(t, x; = Y (t - r, x - £; + J d\ J Y (t - A, x - y; j/)*(A, y; dy; (5.2) 



the functions Zq, Yq were examined in detail in Sect. 4. 

For the functions Q and \l/ we assume the integral equations 



Q(t, x; = M(t, x;0 + JdxJ K(t - A, x; y)Q(X, y; £) dy, (5.3) 



I" 

t 



(t, x- = tf(f , x; + J dxj Kit - A, x; y)#(A, y; dy, (5.4) 
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where 



d 2 

M(t,x;0 = \ K(x) - a i:j ^)]-^-Z (t,x - ^-^) 



dZ Q (t,x-£;€) 
x) + c(x)Z (t,x 



=1 



dxi 



( d 2 

i,j=l ^ 1 3 



+ £ k(x) dYo{t \* Z>® + c(x)Y (t, x - £; 0- 
i=i °^ 



Using the estimates (4.1)-(4.3) we find that 

|M(t, x; 01 < - er n+7 exp{-ap(t, x; 0, 0}- 



(5.5) 



In order to obtain estimates for K, we have to use estimates for Y , different for different 
values of n, and for different domains, and then to roughen the resulting estimates into less 
exact but unified bounds convenient for the Levi method. 



Proposition 4. For any n > 2 

\K(t, x- 1 < Ct^n _ ^\-n+v exp {- (T p(t ) x . 0, £)}, 

TOt/i < 7] < 7. 

Proof. Using Proposition 2 we write the estimate for \x — £| > t a / 2 : 
\K(t,x;0\ < Ct-V+^expi-vpfrx-At)}. 
For I a; — £| < t a l 2 we consider various cases separately. If n = 3, or n > 4, then 
\K(t, x; 1 < Cr^^- 1 |x - £r +2 exp{-ap(t, x; 0,0}- 

Since 

I a: - £| 2 = \x - £\ v \x - t\ 2 - T > <\x- ^t^ 2 ^\ 

(5.7) and (5.8) imply (5.6). Note that transforming (5.7) we change the constant a. 
If n = 2 or n = 4, the initial bound for K is 



(5.6) 



(5.7) 



(5.8) 



|^,x;0l<C7r^-i 



log 



\ x - £1 



t a 



We roughen the estimate (5.9) replacing the factor 
x-^~ n+r > 



{- I ) exp{-<rp(t,a;;0,f)} 



log- 



(5.9) 

+ I with the factor 



t a / 2 



, which results in (5.6). □ 
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Next we have to study the increments 



A x M(t, x- = M(t, x; - M(t, x'; £), 
A x K(t,x;0 = K(t,x;0 ~ K(t,x';0, 



x,x' e R n . 



We can write A x M(t, x; £) = M 1 + M 2 + M 3 where 

m V- A ( ^ 2 Z (t,x-^Q 
M 1 = ^ A xaij (x) — — , 



M 2 = - nij(0 ]A I gM^- r ii« , 



dx\dx'j 



and the term M3 contains lower order derivatives. Since M3 has a weaker singularity and does 
not influence the estimates, we omit its detailed description and consider only M\ and M 2 . 
An estimate for M 1 is given directly: 



I Mi I < d\x - rr'| 7 r a |a; - £|~ n exp{-ap(t,x;0,£)}. 
Using (4.1) we find that 



(5.10) 



A, 



d 2 Z (t,x'-t;0 



dx'fix'j 



< 



d 2 Z 



dx^dx'j 



d 2 Z n 



dxjdxn 



+ 



d 2 Z n 



dx^dx'j 



< Clr™\x - exp{-uap(t, x- 0, t)}\x- x'\ v C\^r a ^ [\x - ^ n exp{-ap(t, x; 0, 0} 

+ K-er ri exp{-ap(t,x / ;0,O}] 1 ^ (5.11) 



where x = x + 9(x' — x), < 9 < 1, < v < 7, and we assumed that min {\x — £\,\x' — £|}^0. 
From (5.10) and (5.11) we find that 

\A x M(t,x;^)\ < C A t~ a \\x - x'p\x - Z\~ n exp{-ap(t,x]0,£)} + \x - x'\ u \x' - ^M(t,x;£;) 

(5.12) 

where 



I l-v 



M(t, x; S) = \x- i\- nv ' v exp{-uap(t, x; 0, 0} [\x - £|~ n exp{-ap(t, x; 0, 0} 

+ |x'-er"exp{- ( Tp(t,x / ;0,O}] 
We will use also the following estimate which is a direct consequence of the definition of A X M: 

\A x M(t, x;0\< Cr a (\x - £P +7 exp{-ap(t, x; 0, 0} 

+ \ x >-Z\-^exp{-ap(t,x'-0,O})- (5.13) 



Let us consider two possible cases. 
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a) Suppose that \x — x'\ > At a / 2 . It follows from (5.13) that 

|A x M(i,a;;OI < C 5 (A)t^^' a [\x - ^\~ n+£ exp{-<7 lP (f, x; 0, 0} 

+\x' -er n+e exp{-(7ip(*,x , ;0,0}] , e > 0, 

whence 

| A x M(i, x; 01 < C 6 r> - x'\^ £ [\x - er n+e exp{-a lP (t, x; 0, 0} 

+|x'-£P +£ exp{-a lP (M';0,0}], e > 0, (5.14) 

b) Let |x — x'| < AW 2 with some A > 0. If \x — x'\ > ^max(|^ — x\, |£ — x'\) with some 
rj > 0, then we obtain again the estimate (5.14). Thus, we have now to consider the most 
complicated subcase, for which 

\x — x'\ < ?7max(|£ — x\, |£ — x'\), r] > 0. 

Here we use the estimate (5.11). 

Without restricting generality, we may assume that |£— x'\ < |£ — x\. Then \x— x'\ <r)\£—x\, 

\x — £| = \x — £ + 6{x' — x)\ > \x — £| — 0r)\x — £| > (1 — r))\x — £|. 

Note that the function 



2 

r \ — 



r^r "exp<J °{—) ?■ ()< n L 

is monotone decreasing. This implies the estimate 

M(t, x;0\< C 7 [\x - t\- n -» exp{-<7ip(t, x; 0, 0} + W - i\' n ~ v exp{-ap(t, x' ■ 0, £)}] ■ 

Setting e = 7 — v and using (5.13), (5.14), we come to the estimate contained in the following 
proposition. 

Proposition 5. If x" is one of the points x,x' for which \x" — £| = vnm.{\x — ^\,\x' — ^\} , then 
\A x M(t,x;0\ < Cr a \x - x'y- £ \x" - £\~ n+£ exp{-o-p{t,x";0,O}, (5.15) 



\A x K(t,x;£)\ < Ct'^x - x'\^- £ \x" - ^\~ n+£ exp{-ap(t,x";0,C)}. (5-16) 

The proof of (5.16) is similar to that of (5.15), since the estimates for Yq have a structure 
similar to those for Zq. □ 

Now we are ready to consider the integral equations (5.3) and (5.4). 

Proposition 6. The integral equations (5.3) and (5.4) have the solutions Q(t,x;£), ty(t,x;£,), 
which are continuous for and satisfy the estimates 

\Q(t,x;0\ < Cr a \x-t;\- n ^exp{-o-p(t,x;0,0}, (5.17) 
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I A x Q(t, x; 1 < Cr a \x - xW' ~ ^ n+£ exp{-<xp(t, x"; 0, 0}, 



(5.18) 



\*(t, x- 1 < Cr 1 \x - £r +7 exp{-ap(t, x; 0, 0}, 



(5.19) 



| A^tf (t, x; 01 < Ct'^x - x'| 7 - £ |x" - £\~ n+£ exp{-ap(t, x"; 0, 0}, (5.20) 

with an arbitrary e G (0,7) and the constants depending on e. If n = 2 or n = 4, £/ien /or any 
p > 

|Q(f , x- 1 < Cr-T | x _ ^j-n+e-/* exp{-ap(t, x; 0, 0}, (5.17') 



\A x Q(t, x- 1 < Cr-" |x - xT V - er" +£ ^ exp{-ap(t, x"; 0, 0}, 

|*(*, x; 1 < Cr 1 -^ \x - £1-^ exp{-ap(t, x; 0, 0}, 

\A x y{t,x;£)\ < Ct- l -^\x - x'\^- £ \x" - er n+e ~ M exp{-o-p(t,x";0,O}- 
i/ere t/ie constants depend on e and p. 



(5.18') 



(5.19') 



(5.20') 



Proof. The equations (5.3) and (5.4) are solved by iterations. The estimates (5.17), (5.17'), 
(5.19), (5.19') follow from Lemma 5. 

Let us prove (5.18). We have A X Q = A X M + J where 



J = J dxj A x K(t-\,x;y)Q(\,y;Ody, 



(5.21) 



t 

\J\ < C|x-x'| 7 " £ J(t- A) _1 A _Q dA J \x-y\- n+£ \y-£\- n+ ^exp{-ap(t,x;\,y) 
B" 

- <rp(A, y- 0,Z)}dy + J \x' - y\- n+£ \y - exp{-ap(t, x'; A, y) - ap(X, y- 0, 0} dy 

R" 

t r 

< C\x-x'\^- £ J(t- A)- 1+ ^ £ A (1 - Q) - 1 rfA y |x-y|-^%-f|- n ^exp{-(7ip(t,a;;A,j/) 
un« 

- o- lP (A, y; 0,Z)}dy+J \x' - y\~ n+ i \y - exp{-<7 lP (t, x'; A, y) - a lP (X, y; 0, 0} dy 

By Lemma 3, the integral (5.21) is estimated as follows: 

|J| < C\x - x'\i- e r a+ ^ [|x-er n+7+f exp{-o- 2 p(t,x;0,O} 

+ |x / -er" +7+f exp{-a 2 p(t,x;0,O}] 
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Together with (5.15), this implies (5.18) with an appropriate a > 0. 
The proof of the inequalities (5.20), (5.20') is similar. □ 

5.2. Heat potentials. In order to verify that the functions (5.1) and (5.2) indeed solve the 
Cauchy problem via the representation (1.4), we have to study the heat potential 

t 

W{t, x) = J d\J Y (t - A, x - y; y)f(X, y) dy (5.22) 

R« 

for two situations. Firstly, we must consider the case of a bounded, locally Holder continuous 
(in x) function /, jointly continuous in (t,x) G [0,T] x W n . This is, of course, dictated by 
the formula (1.4) containing a term which solves the inhomogeneous equation with the zero 
initial condition. Note that Yq is the main singular part of the function Y appearing in (1.4). 
Secondly, heat potentials appear also in (5.1) and (5.2), where the role of / is played by the 
functions Q and \l/ which are much more singular (see Proposition 6). 

Let us consider the potential (5.22) with a bounded, locally Holder continuous /. The 
existence of the integral (5.22) and the possibility to find its first order derivatives in x by 
differentiating under the integral's symbol follow directly from the estimates (4.5), (4.7), (4.9), 
and (4.11). For studying the second order derivatives we consider the function 

t-h 

W(t, x) = J dxj Y (t — X,x — y; y)f(X, y) dy, h > 0. 

R" 

As t — X > h, the function Y Q and its first and second derivatives have integrable singularities 
at x — y. Therefore we may differentiate under the integral's symbol, so that 

d 2 W h (t,x) _ 

— J\ + J2 



dxidxj 



where 

t-h 



Ji = JdX / 9y ° (t ^," y;y) [/(A,y)-/(A,x)]dy, (5.23) 

R n 



We decompose further, J\ = + j[ 2 \ where corresponds to integration over the 
set = \y G W 1 : \x — y\ < (t — A) a / 2 } while for j[ 2 ^ the domain of integration is = 
[y G W 1 : \x — y\ > (t — A)"/ 2 }. Suppose that, for example, n > 4; other cases are treated 
similarly. If y G n^, then 



d 2 Y (t- X,x-y;y) 



dxidxj 



< C(t- xy a - 1 \x- y \ 
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so that the part of the integral in (5.23) is majorized by 



t-h 



li 1] =C J(t-X)- a ^d\ J \x - y\- n+2+ ^ dy, 



o n(i) 
and the change of variables y — (t — X)~ a l 2 (y — x) gives 



t-h 



I?<C j(t-Xr 1+ ^d\ J \y\ 

|y|<i 



-n+2+7 



dy. 



(5.25) 



Here 7 > is the Holder exponent of the function f(X,x) in x. 
For y e Il {2) we use the inequality 



d 2 Y (t- \,x-y;y) 



< C(t- A)-^" 1 exp^ -a 



\x - y\ 
(t - A) Q A 



a consequence of the general estimate (4.11). A similar argument shows that the n^ 2 ) part of 
the integral in J 1 is majorized by 



t-h 



C J (t - A)- 1+ ^ dX J |y| 7 exp{-<j|y|^} dy. 



|£|>1 

We see the existence of the limit of J\ as h — > 0. 
Let us consider J 2 . Using (4.14) we have 



t-h 



J2 



f(X,x) dX 



d 2 Y (t - X, x - y; y) d 2 Y (t - X, x - y; 5) 



dxidxj 



dxidxj 



&=x. 



dy. (5.26) 



By (4.12), the expression in brackets is majorized by 



C\x - y\- n+2+ \t - A)-"" 1 exp{-ap(t, x; A, y)}. 
The same change of variables (used above) shows that the integrand in (5.26) is majorized by 



(t-X) 



-!+l?|y|-«+ 2 +7 



cxp j— a\y\ 2 -« I . 



This means the existence of the limit of J 2 as h — > 0. Simultaneously we have proved the 
formula 



d 2 W(t,x) 
dxjdxj 



= / dX 



d 2 Y (t- X,x-y;y) 



dxidxj 



[f(X,y)-f(X,x)]dy 
t 

d 2 Y {t- X,x-y;y) 



+ 



J f(X,x) dX J 



dxidxj 



dy. (5.27) 
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Now we have to carry out a similar procedure for the potential (5.22) with f(X,y) = 
Q(X,y;£) with a fixed £ (the case of ^(A,y;r, £) is quite similar). We use the estimates (5.17) 
and (5.18) for the function Q and its increment, and also the following roughened estimate: 



d 2 Y (t-X,x-y;O 



dxidxj 



<C(t — X) \x — y\ n exp{-ap(t,x; X, y)}. 



(5.28) 



In fact we have to prove convergence of the integrals 



t 



d 2 Y (t- X,x-y;y) 



R" 



dxidxj 



\Q(\,x;Z)-Q(\,y;0\dy, 



r' 2 = J \Q(X,x;0\dX J 



d 2 Y (t - X, x - y- y) d 2 Y (t - X, x - y; 5) 



dxidxj 



dxidxj 



S=x 



dy. 



Let a > be a small positive constant. Changing a we can rewrite (5.28) as 
d 2 Y (t-X,x-y;O 



dxidxj 



< C(t- X)- 1+ ~\x - y\~ n - a exp{-ap(t,x; X,y)}. 



Now we get by Lemma 3 and Lemma 4 that for < e < 7 



z 

r[<C J{t-X) 



-i+- 



fx-«dx 



I 



x-y 



—n+i—s—a 



n+£ exp{-ap(t,x; X,y) 



-ap(X, x- O,0}dy + J\x- y\- n+ ^ £ ' a \y - £P +£ exp{-ap(t, x; A, y) - ap(X, y; 0, £)} dy 

t 

= C\x - £\- n+£ exp{-ap(t, x- 0, 0} J (t - Xy 1+2 ?X- a dX 



x / \x — y\ 



-n+y—e—a 



exp{-ap(t,x;X,y)}dy 



1 

+ C j(t-Xy 1+ ^X- a dX I \x-y\- n+ ^- a \y-C\- n+e ex P {-ap(t,x;X,y)-ap(X,y;0,O}dy 



< C\x - i\~ n+£ exp{-crp(A,x;0,£)} / (t - X^^^^X^dX 



+ Ct- a+a -?\x - Z\- n ^- a exp{-ap(t,x-0,O}- 



As before, the letters a, C meant various positive constants; this abuse of notation will be 
convenient in the sequel too. 



30 



In order to obtain a bound for I' 2 , we use the estimate 
d 2 Y (t - X,x - y;y) d 2 Y (t - X,x-y;S) 



dxidxj 



dxidxj 



S=x 



< C\x - yy- n {t - X)- 1 exp{-ap(t, x; X, y)} < C\x - y\~ n+ ~i{t - A)~ 1+ ^ exp{-ap(t, x; X, y)} 
(with a different a). Then by Lemma 4 

I' 2 < C\x - £|- n+7 exp{-ap(t, x; 0, £)} 
t 

x J(t- Xy 1+ ^X- a dX J \x-y\- n+ % exp{-ap(t,x;X,y)} 

I" 

t 

< C\x - er ri+7 exp{-ap(t, x; 0, £)} J (t - Xy l+ ^X~ a dX 

o 

< C t^- a \x - £r +7 exp{- < rp(t, x; 0, 0}- 
Thus, we have proved (5.27) for f(X,y) = Q(X,y;£). 

5.3. TTie fractional derivative. Our next task is to study the fractional derivative of the 
heat potential (5.22), again for the above two cases. 

Suppose that f(t,x) is bounded, locally Holder continuous in x, and jointly continuous in 
(t,x). We have O^VF = J^t> where v = I^ a W is the Riemann-Liouville fractional integral, 
that is 

t e 

v(t, x) = 1 Jit- 9)- a dejdxj Y (9 - X, x - y- y)f(X, y) dy 



o 

t t 



j^—^ jdxj(t- ey a dejY (9-x,x- y - y )f(x, y ) 



dy 



o A 



dX J Z (t - X,x - y;y)f(X,y)dy, 



R n 



by the definition of Y . 

Let h be a small positive number, 



t-h 



v h (t,x) = J dX J Z (t - X,x -y;y)f(X,y)dy. 

R n 

Then v h — * v pointwise as h — > 0, 



= «?>(*,*) +«?>(«,*) 
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where 



f(Ky) dy. 



1" 



We have 



v^\t,x) = J [Z (h,x -y;y) - Z (h,x - y;x)]f(t - h,y) dy 

i" 

+ J Z (h,x-y;x)[f(t-h,y)-f(t-h,x)]dy + f(t-h,x), (5.29) 

due to (4.13). By Proposition 3, the absolute value of the first integral in the right-hand side 
of (5.29) does not exceed 



CK 



a J |x-y|- n+2+7 exp|-(7 (h- a/2 \x - y\) 2 - a } dy = const • -> 0, 



as h — > 0. Similarly, using (4.1) we show that the second integral in (5.29) tends to zero as 
h -> 0. Thus 

v£\t,x) — /(t,x), (5.30) 

as h — > 0. 

Turning to t> ^ we write 

<9Z (t -\x-y;y) 



v£\t,x)= JdxJ 



dt 



-[f(X,y)-f(X,x)]dy 



R" 



t-h 



, . 8Z (t- X,x-y;y) 
I f(X,x)dX I — dy. 



It follows from (4.15) that 



t-h 



I dX I 



dZ (t -X,x-y;y) 
dt 



[f(X,y)-f(X,x)]dy 



<c(t-xy 1+ ^. 



Together with (4.16) this shows that 



R n 



dZpjt - X,x - yyy) 
dt 



[f(X,y)-f(X,x)]dy 



+ 



J f(X,x)dX J 



dZpjt -X,x-y;y) 
dt 



dy. (5.31) 
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dZ 

By the definition of Yq, — = Dj Yq, and it follows from (5.30) and (5.31) that the 

L/ v 

fractional derivative D^jy exists and can be represented as follows: 

t 

(d?V) (t,x) = f(t,x) + jdxj dZo{t ~ X m X ~ V] V) [f{\ V) - f{\ x)\ dy 



MP 



+ jf{\x)d\j dZ ° (t -^- mV) dy. (5.32) 



The case, in which /(A, y) = Q(X, y; £), is considered similarly, on the basis of the estimates 
(4.15'), (4.16), (5.17), (5.18). 

5.4. The initial condition. The above study of heat potentials, together with the investiga- 
tion of the integral f Z (t, x — £) M o(0 d£ performed in jT3| shows that our construction of the 

R™ 

kernels Z, Y indeed gives, via the formula (1.4), a solution of the equation (1.1). It remains to 
verify the initial condition (1.3). 
By our construction, 

Z(t, x; i) = Z {t, x-C,0 + V z {t, x- 0, 

where 

t 

V z (t, x- = I dX Y (t - X, x - y; y)Q{X, y; dy 



o 



(see (5.1)). We have 
J Z (t,x - £,£)«q(0 d£ = J Z (t,x-£,x)uQ(£)d£ 

R" R" 

+ J [Z {t, x - e, - Z (t, x-£, x))u {0 

The first summand tends to Uo(x), as t — > [T^]. By Proposition 3, the second summand is 
majorized by 

Ct' a J |s-^|- n+2+7 exp{-a (r a/2 \x - d£ 

MP 

= CV a J \r]\- n+2+ ~< exp {-o-\7]\^} dq — ► 0, 

i™ 

as t — > 0. 
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Let us find an estimate for Vz- Suppose, for example, that n = 3 (other cases are treated 
similarly). Then by (4.7) and (5.17) 



\Vz(t,x;£)\<C j{t-\)-^ l \~ a d\ j \y-^\- 3 ^ex V {-a[p(t,x;\,y)+p(\,y;0,O}}dy. 

R3 

Taking < 70 < 7 we can rewrite this as 

\Vz(t,x;t)\<C [(t-Xr^X^dX [ ( J g ~ y J /9 V +7 ° \x - yj- 1 -^ - ^ 







(t- X)«/\ 
x exp{-a[p(t, x; A, y) + p(X, y; 0, f )]} dy 



<Ci jit-Xy'^X^dX J \x-y\- 1 -^\y-Z\- 3 ^eM-°Mt,x ] X,y)+p(X,y-,0,O}}dy, 

R3 

and now Lemma 3 implies the estimate (iii) from the main Theorem. 

Other estimates given in the formulation of our main Theorem are proved in a quite similar 
way (note that we roughen some estimates containing logarithmic terms replacing the logarithm 
by the power function with an arbitrarily small exponent). Returning now to the verification 
of the initial condition we obtain that 

J V z (t,x;0M0d£ <Ct^- a j ft™- 1 exp {-a (t-*' 2 ]^} 

= Ct°^ J l^l 7 " 70-1 exp j-o"M^j dr] — ► 0, 

as t — > 0. Thus, the initial condition (1.3) has been verified. 

5.5. The Holder continuity in t. Let us prove the assertion b) of the main theorem. Suppose 
that < v < I — a, a + X < v. We have to show that for any fixed 16K" t u u(t, x) is Holder 
continuous in t with the exponent a + A. It is sufficient to prove the Holder continuity near the 
origin t — 0. Let us consider, for example, the case n = 3; all other cases are treated similarly. 
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Let us write u(t, x) = u±(t, x) + u 2 (t, x) + u 3 (t, x) + w 4 (t, x) where 

u 1 (t,x) = J Z (t,x - C;0 M o(0 d£, 

u 2 {t,x) = I V z (t,x;£)u (O d £, 



z 

u 3 (t,x) = J dxj Ko(i-A,x-£;£)/(A,Od£, 





t 



u 4 (t,x) = J dXj V Y (t-\,x;0f(K0d£. 

R3 

The Holder continuity of t v U\ (t, x) was proved in [T3| (with the use of fractional calculus 
and the asymptotics of H-functions) for the case of constant coefficients. The proof for our 
case, in which there is also a dependence on the parameter £, is identical to [13 . 

As we saw above, \u 2 (t,x)\ < Ct a ^/ 2 , so that \t u u 2 (t,x)\ < Ct^ +U < Ct a+X for small values 
of t, since a + X < v. 

Next, by (4.7) 

t , ^_>. t 



\u 3 (t,x)\<C Jit-Xy^dX J expj-<r( (t _ A) J /2 1 M£ < C x J (t - X)^ 1 dX = C 2 t c 

R3 ^ ^0 

so that |t I/ M 3 (t,x)| < C 2 t a+V < C 2 t a+X for small values oft. 

Finally, from the estimate (iii) of the main Theorem we find that 



\t v u 4 (t,x)\ < Ct^ +a+u < Ct 



a+X 



for small values of t, as desired. 



6 THE ONE-DIMENSIONAL CASE 

In the case n = 1 the parametrix and the fundamental solution have no singularity in the 
spatial variable. This simplifies the situation greatly and makes it similar to the theory of 
conventional parabolic equations of an arbitrary order. Note however that even this case is 
more complicated than the classical study of a second order parabolic differential equation. In 
particular, the estimates of iterated kernels in the Levi method should still be performed in two 
stages (see Lemma 6). 

Here we give the main estimates for this case; their proofs, simplified versions of those given 
the preceding section, are omitted. The scheme of the Levi method and the main notations are 
as above; see (5.1)- (5. 4). 
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Proposition 7. If n — 1, then 

\M(t,x;0\ < Cr^ a / 2 exp{-ap(t,x-0,O}, 

\Q(t,x;0\ < Cr( 3 -^/ 2 exp{-ap(t,x;0,O}, 

The estimates for differences of the above kernels are collected in the next proposition. 
Proposition 8. If n = 1, then 

\A x M{t,x;0\ < C\x - i'|n-^J exp{-ap(t,x";0,£)}, 
\A x K(t, x;0\< C\x - x '|7-e r i-(i- e )a/2 expj-ap^, x "; 0, 0}, 

\A x Q(t, x;£)\< C\x - x'\^ £ t- {3 ~ £)a/2 exp{-ap(t, x"; 0, 0}, 
\A x q{t,x;0\ < C\x - x'l^-H- 1 ^ 1 -^ 2 exp{-ap(t,x";0,O}, 

< e < 7. 



7 NONNEGATIVITY 

Let us prove that the functions Z and Y are nonnegative. It is sufficient to show the non- 
negativity of the solution (1.4) of the problem (1.1), (1.3) for arbitrary nonnegative functions 
uo e tf °°(R n ), /eC o °°([0,T] xJET). 

As we already know, the solution u(t, x) is bounded (due to our estimates of Z and Y), say 
\u(t, x)\ < M, and belongs to if" +A [0, T], < A < 1 — a, a + A < u, for each x. This regularity 
in t makes it possible to represent the regularized fractional derivative B^ U via the Marchaud 
fractional derivative: 

t ] u) {t,x) = 1 [r a u(t,x) - r a u {x)] + a HmJ £ {t,x) 
/ 1 (1 — a) 1 (1 — a) e^+o 

where ip e (t, x) = for < t < e, 

t-e 

, , . f u(t,x) — u(t,x) , _ 



and the limit exists for all x G M n , t G [0, T] (see Sect. 13 in [2T]). 

Denote v(t,x) = — 7^^\ "where is the Mittag-Leffler function 0, /3 > will be chosen 

later. Using the fact that B[ a) E a ((3t a ) = f3E a (f3t a ) we obtain that 

B { t a) u(t,x) = (3E a {(3t a )v{t,x) - — \ ru (x) + (L p v) (t,x) 

L [1 — a) 
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where 



(L p v) (t,x) 



r(i-a) 



t-e 



v(t,x)r a + a lim / E a (Pr a ) V ^ X) ^^- dr 



Now the equation (1.1) takes the form 

1 



where 



E a (pt° 
g(t,x) 



(Lpv) (t, x) - (B - @)v(t, x) = g(t, x) 



(7.1) 



E a (/3t<*) 



r(i-a 



-u (x) + f{t,x) 



Let Co = sup |c(x)|, d > 0, (3 — cq + d. Consider the function 

F R (t,x) = — (\x\ 2 + fit a + 1) , R,fi>0. 
Since F R (t,x) and E a (f3t a ) are monotone increasing in t, we have 

(W (t, X ) > W+ F R ) (t,x) > ^D° 0+ (t°) = jW±g> . 

Here -Dq + is the Riemann-Liouville fractional derivative (see [21]), and in estimating LpF R we 
used the Marchaud form of Dq + . 
On the other hand, 



n 2M n 2M 

(B - /3)F R (t,x) = 22au(x)-jp + 22 b j( x )-^ x j + i c ( x ) ~ co\F R (t,x) - dF R (t,x) 



3=1 



R 2 



M 



<^(Ci + C 2 \x\ - d\x\ 2 - d^t a - d) 
R 2 



so that 
1 



M 



j, w (W(u)-M%«)> ? 



A*r(l + q) 



Ci - C 2 |x| + d|x| 2 + + d 



for all s G K n , t G [0, T], if the number fi is large enough. 

Let G(t,x) = v(t,x) + F R (t,x). Since g(t,x) > 0, it follows from (7.1) and (7.2) that 



(L p G) (t, x)-{B- p)G(t, x) > 0. 



> 
(7.2) 

(7.3) 



If \x\ = R, then 



G(t, x) = x ) + M + + — > 0. 



R 2 R 2 

Also G(0,x) = u (x) + F R (0,x) > 0. Then G(t,x) > for all (t,x) with t G [0,T], |ar| < R. 
Indeed, otherwise the function G possesses a point (t°,x°) of its global minimum on the set 
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{0 < t < T, \x\ < R}, such that G(f,x°) < 0. At this point (B - f3)G(t°,x°) > 10J, and by 
(7.3) (LpG) (t°,x°) > 0, which contradicts the definition of hp. 
Thus, if \x\ < R, then 

u(t, x) > (\x\ 2 + fit a + 1) E a (f3t a ). 
Since R is arbitrary, this means that u(t,x) > 0. 
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